Interfaces of modulated phases 
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Numerically minimizing a continuous free-energy functional which yields several modulated 
phases, we obtain the order-parameter profiles and interfacial free energies of symmetric and non- 
symmetric tilt boundaries within the lamellar phase, and of interfaces between coexisting lamellar, 
hexagonal, and disordered phases. Our findings agree well with chevron, omega, and T-junction 
tilt-boundary morphologies observed in diblock copolymers and magnetic garnet films. 



Modulated phases are found in a surprisingly diverse 
set of physical and chemical systems, including super- 
conductors, thin- film magnetic garnets and ferrofluids, 
Langmuir monolayers, and diblock copolymers |J. Such 
phases are characterized by periodic spatial variations of 
the pertinent order parameter in the form of lamellae, 
cylinders, or cubic arrangements of spheres or interwo- 
ven sheets. This self-organization results from competing 
interactions: a short-ranged molecular one favoring a ho- 
mogeneous state, and a long-ranged contribution, which 
can have magnetic, electric, or elastic origin, favoring 
domains. Because of the modulation, interfaces between 
different phases or grain boundaries within a single phase 
are most interesting, yet they have received much less at- 
tention than those occurring in solids. Recently, experi- 
mental studies of grain boundaries within lamellar phases 
of diblock copolymer have been carried out (||| which 
illustrate the rich interfacial behavior exhibited by such 
systems. In particular, three morphologies of tilt bound- 
aries[2b], denoted chevron, omega, and T-junction, were 
observed. Such interfaces are more difficult to describe 
theoretically than those between uniform phases, which 
have been the subject of classic work jij. In this article 
we study not only tilt grain boundaries within lamellar 
phases, but also interfaces between coexisting modulated 
phases of different symmetry. 

The dimensionless free-energy functional we use, 
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includes an enthalpic term (proportional to the interac- 
tion parameter x)i favoring an ordered state in which 
\4>\ is non-zero, an entropy of mixing preferring a dis- 
ordered state, (j> — 0, and confining |</>| to be less than 
unity, and derivatives of the order parameter. The or- 
dered state occurs with a modulation of dominant wave 
vector q* = 1/V2 because of the competition between 
the negative gradient square term (favoring domains at 



large length scales) and the positive Laplacian square 
(preferring a homogeneous state at small length scales). 
Such a free energy functional has been used to describe 
the bulk phases of magnetic layers Q, Langmuir films 
amphiphilic systems, jjj, and diblock copolymers [||, 
and the effects of surfaces on isotropic M and hexagonal 



phases 10 of the latter. 



In our numerical studies, we determine the minimum 
of (1) directly using a conjugate-gradient method, which 
is more convenient than solving the corresponding Euler- 
Lagrange equation. We discretize on a square lattice, 
approximate the derivatives by nearest-neighbor differ- 
ences, and employ a mesh size sufficiently small that 
discretization effects are negligible within the numerical 
precision. A typical grid contains 40,000 points. Fig.l 
shows the resulting two-dimensional bulk phase diagram 
as a function of x and a) the average order parameter 
<f> = ((f)) v and b) the chemical potential /i. It is in good 
agreement with previous calculations based on single- 
mode approximations j5|,[| . In addition to the disordered 
(D) phase, there is a lamellar (L) and two hexagonal (H) 
phases, which all join at the critical point at \c = 3/4, 
fj, = $ = 0. (In three dimensions, one expects additional 
cubic phases, which are not studied here.) For larger val- 
ues of Xi f ne H an d L phases each terminate at a triple 
point, whereas experimentally one sees that these mod- 
ulated phases exist even for very large values of x- This 
unphysical part of the phase diagram is due to a break- 
down of the gradient expansion in (1). 

We first present results for grain boundaries in lamellar 
phases, and begin with the asymmetric tilt grain bound- 
ary (GB) between two perpendicular L phases. This is 
the T-junction of Ref . [2b] , for which the layer continuity 
between the two adjoining phases is disrupted. Figs. 2a- 
c show order-parameter profiles <j) f° r different x. an d 
fi = (symmetric stripes), while Figs.2d-f show profiles 
for varying chemical potential /1 and x — 1 (asymmet- 
ric stripes). Fig. 2b clearly shows the enlarged endcaps 
noted in experiment and the series 2a-c predicts 

that they become less pronounced with increasing x- The 
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GB interfacial energy job scales as job ~ (x ~ XcY* 
with yu % = 3/2, Fig.2g, in accord with mean-field pre- 
dictions jllj . To demonstrate that the interfacial energy 
7gb is indeed quite small, we show with a broken line 
the difference in bulk free energies between the disor- 
dered and lamellar phases multiplied by one wavelength 
of the latter. This is roughly the cost per unit area of 
disordering such a width of lamellar phase. Clearly the 
actual grain boundary bridges the two grains in a man- 
ner much less expensive than the insertion of a region 
of disorder. In Fig.2h we plot the interfacial energies 
of two distinct GB structures, which are degenerate at 
(i = 0, being related by order-parameter reflection sym- 
metry — > — 0. For /i =/= this symmetry is broken and 
the two structures correspond to distinct local free en- 
ergy minima, one metastable, the other globally stable. 
The free-energy barriers between such interfacial struc- 
tures are responsible for slow interface motion and thus 
long healing times in multi-grain lamellar samples |Q ] . 

In Figs.3a-c we show symmetric tilt-boundary (TB) 
configurations for a fixed angle 8 = 90° between the 
layer normals as a function of x- Here the layer continu- 
ity is maintained across the boundary, and the chevron 
morphology is quite evident. In Figs.3d-f, the tilt an- 
gle is progressively increased at fixed X-, an d one clearly 
sees the change from the chevron to the omega structure. 
The omega shape of the layers at the boundary results 
from frustration due to an imposed local lamellar wave- 
length much larger than the equilibrium value. Figs. 3d- 
f resemble micrographs of undulating lamellar patterns 
in garnet films 113], and the similarity of Fig.3f to the 
micrographs [2b] of diblock copolymer TB's is striking. 
The scaling of the TB energy jtb is again described 
by an exponent /u* = 3/2 (Fig.3g). Close to critical- 
ity, one finds pronounced reconstruction in terms of a 
square-like modulation, Fig. 3a. For small 8, Jtb{0) ~ 8 3 , 
Fig.3h, in accord with the bending behavior of elastic 
sheets; the 8 — > 180° limit is expected to be linear, 
1tb{8) ~ 180° — 8, in accord with a description in terms 
of decoupled dislocations of finite creation energy. 

We now turn to interfaces between thermodynamically 
distinct phases. Interfacial structures between coexisting 
lamellar and disordered phases for three different values 
of the angle ■& between the lamellae and the interface 
are shown in Figs.4a-c for x = 1 (where the phases and, 
therefore, also the interfacial structure are metastable), 
and in 4d-f for \ — 1-5 (where they are stable). The 
metastable L-D boundary is shown as a dashed line in 
Fig. lb. In Figs. 4a and d, for # = 0°, there is a relaxation 
of the outermost layers as the interface is approached, 
leading to a small increase in the wavelength, as in solids. 

Before presenting the interfacial free energies, we de- 
scribe our method of obtaining them, because the calcula- 
tion of such free energies between coexisting phases which 
can be modulated is non-trivial |L4j . We calculate the to- 
tal free energy, F 1 , of phase J in a box employing periodic 



boundary conditions parallel to the left and right faces, 
and reflecting (Dirichlet) boundary conditions on those 
faces themselves. We adjust the length of the box so 
that the free energy is minimized. This occurs when the 
length of the box is some integer number of wavelengths 
of the periodic structure of phase /. The volume of the 
box is V . By this means, there is no surface contribution 
to the total free energy, so that the bulk free energy is 
obtained directly; // = F 1 (V')/V . In a similar way we 
obtain f^ 1 = F H (V")/V". For the system in I, II coex- 
istence, we calculate the total free energy in a box large 
enough that the order parameters attain their bulk values 
on the left and right faces at which reflecting boundary 
conditions are employed, and we vary the length of the 
box to minimize the total free energy. The volume of the 
box is V. Again there are no surface contributions so 
that we obtain F^V.A) = V(f£ + / 6 ")/2 + A luI . As 
the bulk free energies are known, the desired interfacial 
energy follows directly. 

The interfacial energy jld is shown in Fig.4g. As x is 
decreased, the coexistence between L and D is preempted 
by the hexagonal phase (see Fig. lb). This is manifest in 
the behavior of 7ld for d = 90° (broken line in Fig.4g) 
which becomes negative. This value of ^ld is obtained by 
assuming an interfacial reconstruction locally resembling 
the energetically preferred hexagonal phase (see Fig. 4c). 
The value of jld for $ = 0° (solid line in Fig.4g) re- 
mains positive down to \ — 0-82 and vanishes there with 
the classical tricritical exponent /i* = 2 p}]. Strong 
anisotropy of as function of d is observed (Fig.4h), 
from which the shape of lamellar droplets in an isotropic 
phase can be qualitatively inferred: For \ > 1-34 the 
drops are elongated parallel to the lamellae; for \ < 1-34 
the L and D phases are metastable at coexistence and 
the drops are elongated perpendicular to the lamellae. 

Hexagonal-disordered interfacial structures are depic- 
ted in Figs.5a-c; the corresponding interfacial energy 
7hd, shown in Fig.5g, scales again with /i* = 2. Struc- 
tures for the hexagonal-lamellar interface are plotted in 
Figs.5d-f, where we chose the mutual orientation in ac- 
cord with the epitaxial relationship found for amphiphilic 
systems [^6|. The corrugation of the lamellae near the 
interface, particularly evident in Fig.5d, resembles that 
seen in experiments on diblock copolymer blends p7| . 
The interfacial free energy of this interface, jhl, is plot- 
ted in Fig.5h and scales with the classical critical expo- 
nent //* = 3/2. This reflects the fact that, in contrast 
to the L-D and H-D interfaces, the H and L phases are 
locked into a fixed relative position with respect to trans- 
lations perpendicular to the H-L interface. 

At the triple point between disordered, lamellar, and 
hexagonal phases, x — 1-34, we find 7ld > Ihl + Ihd- 
It follows that the L-D interface that we have calculated 
is not the thermodynamically stable one. Within our 
model, the L-D interface is therefore wetted by the hexag- 
onal phase at the triple point. As the occurrence of such 
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points between these three phases is not uncommon ex- 
perimentally Jl8| , it would be interesting to determine 
whether this wetting does indeed occur. 

Density fluctuations will change the interfacial criti- 
cal exponent from /i, = 3/2 to /i* = 1.26, but leave 
the classical tricritical exponent /x* = 2 intact jllj , while 
fluctuations of the direction of the modulation normals 
will eliminate the critical point p9[ . By introducing uni- 
axiality, these latter fluctuations can be suppressed. 

In summary, we have employed a simple Ginzburg- 
Landau free-energy functional to calculate profiles and 
free energies of several interfaces of modulated phases. 
Qualitative agreement with experiment is very good. The 
observed chevron and omega morphologies at lamellar 
grain boundaries emerge naturally, as do the expanded 
endcaps characteristic of the T-junction. We also cal- 
culated profiles and free energies of interfaces between 
disordered and modulated phases, and between modu- 
lated phases of different symmetry. In all but the sim- 
plest cases, there is significant reconstruction which leads 
to low interfacial energies. An extreme example of such 
reconstruction is the lamellar-disordered interface, which 
we find to be wetted by the hexagonal phase at an L-H-D 
triple point. 
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FIG. 1. Two-dimensional bulk phase diagram, showing 
disordered (D), lamellar (L), and hexagonal (H) phases, as 
a function of the interaction strength y an d a ) the average 
order parameter $ and b) the chemical potential fi. Broken 
lines in a) denote triple lines, broken lines in b) denote the 
(metastable) L-D transitions which exhibit tricritical points 
(denoted by filled circles). 

FIG. 2. Grain boundary (GB) between two perpendicular 
lamellar phases, a-c) Contour plots of the order parameter 
profiles for [i = and y = 0.78, 1, 1.5; throughout the ar- 
ticle the order parameter range [—1,1] is represented by 20 
grayscales; d-f) Profiles for y = 1 and ^ = -0.02, 0.02, 0.06. 
g) GB interfacial energy ycB for jj, — as a function of y, 
showing asymptotic scaling with an exponent /i» =3/2 (in- 
set). The dashed line gives the surface free energy density 
of one lamellar layer for comparison, h) "/gb for y = 1 as a 
function of [i for two GB configurations which are degenerate 
at jj, = 0, but distinct otherwise. 

FIG. 3. Tilt boundary (TB) between two lamellar phases, 
a-c) Profiles for fj, = 0, a fixed angle 6 = 90° between the 
layer normals and y = 0.78, 1, 1.5; d-f) Profiles for /i = 0, 
X = 1, and 61 = 28.08°, 53.14°, and 126.86°. g) TB interfacial 
energy 7ts for /i = as a function of y, showing asymptotic 
scaling with an exponent ^t, —3/2 (inset), h) 7tb for y = 1 
as a function of 9, showing a 7tb ~ 9 3 behavior for small 8 
(broken line) and linear behavior for 9 — > 180° . 

FIG. 4. Lamellar-disordered interface, a-c) Profiles for 
X = 1 and different angles # = 0°, 45°, and 90° between 
the lamellae and the interface; d-f) profiles for y = 1.5. g) 
Interfacial energy jld for "& — 90° (broken line) and # = 0° 
(solid line), the latter scaling with (i, = 2 on approach to the 
metastable tricritical point (inset), h) jld as a function of 
the angle for y = 1 (broken line) and for y = 1.5 (solid 
line). 

FIG. 5. a-c) Hexagonal-disordered interfacial profiles for 
y = 0.78, 0.9, and 1.2; d-f) Hexagonal-lamellar profiles for 
y — 0.78, 0.9, and 1.2. g) Interfacial energy ^hd along HD 
coexistence as a function of y, scaling with /i* = 2 (inset); h) 
7hl as a function of y, scaling with [i* = 3/2 (inset). 
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